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Abstract 

A method is described which, under the assumption of SU(3) symmetry, allows 
one to determine the angle 7 = Arg(V* b ) of the unitarity triangle from time- 
independent measurements of the branching ratios for the rare two-body decays 
B + — * tt°K + and B~ — ► ir K~, as well as of the CP-averaged branching ratios 
for the decays i? ± — > tt^K and — * Tr^-zr , all of which are of order 10^ 5 . The 
effects of electroweak penguin operators are included in a model-independent way, 
and SU(3)-breaking corrections are accounted for in the factorization approxima- 
tion. 
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The study of CP violation in the weak decays of B mesons will provide important tests 
of the flavor sector of the Standard Model, which predicts that all CP violation results 
from the presence of a single complex phase in the quark mixing matrix. The precise 
determination of the sides and angles of the unitarity triangle, which is a graphical 
representation of the unitarity relation V* b V u d + V* b V c d + V t * b V t d = 0, plays a central 
role in this program I0J. Whereas the angle (3 = —Arg(V t d) will be accessible at the 
first-generation B factories through the measurement of CP violation in the decay B 
J/ipKs, the angle 7 = Aig{V* b ) is harder to determine. The sum (P+j) can be extracted 
in a theoretically clean way from measurements of CP violation in the decays B — > tttt 
(or in the related decays B — > irp and pp), but because of experimental difficulties such 
as the detection of the mode B — > 7r°7r° this will be a long-term objective. A method to 
determine 7 proposed by Gronau and Wyler uses rate measurements for six B — > DK 
decay modes ||, some of which require the reconstruction of the neutral charm- meson CP 
eigenstate D° + . A variant of this approach using B — > DK* decays has been discussed by 
Dunietz |3j . Unfortunately, these methods rely either on measurements of some processes 
with very small branching ratios, posing experimental HI] and theoretical 0] challenges, 
or on measurements requiring considerable precision (see, e.g., Refs. || [7j and references 
therein). 

In view of these difficulties, approximate methods to determine the angle 7 have 
received a lot of attention. The simplest of these methods was proposed by Gronau, 
Rosner and London (GRL), who suggested a triangle construction involving the ampli- 
tudes for the decays B + — > tt + K°, tt°K + , and tt + tt°, as well as for the corresponding 
CP-conjugated decays ||. Besides a plausible dynamical assumption this method relies 
on SU(3) flavor symmetry in relating B — > tttt with B — > ttK decays. Later, it was ar- 
gued that the GRL method is spoiled by electroweak penguin contributions, which have 
an important impact in B — > ttK decays and upset the naive SU(3) triangle construc- 
tions H [Tc|l . More sophisticated methods based on quadrangle constructions involving 
other decay modes such as B° s — ► TT°r] p| or B + -»• rj®K + P, pi were invented to 



circumvent this problem. There have also been proposals for deriving bounds on 7 us- 
ing CP-averaged rate measurements in B — ► ttK decays JT3], [II, [T3|, [T(| [T7], and 



for combining these measurements with those of rate asymmetries and other decays like 



B — > KK to obtain further information ITSI, [20 



In the present note we propose a variant of the original GRL method, which based 
on the findings of our previous work flTj| includes the potentially dangerous electroweak 
penguin contributions in a model-independent way using Fierz identities and SU(3) sym- 
metry. We thus obtain an approximate method for learning cos 7 that is conceptually as 
simple and uses the same experimental input and theoretical assumptions as the GRL 
method, though the actual triangle constructions are somewhat more complicated. The 
main advantage of our approach is that it is based on rare two-body decays that are 
relatively easy to access experimentally, and that have larger branching ratios than the 
decays needed for all other methods of measuring 7. Although the accuracy of this 
extraction may ultimately be limited by theoretical uncertainties, even an approximate 
value for cos 7 will be very useful, if only to help eliminating discrete ambiguities inherent 
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in other determinations . 

The basis of our method is the amplitude relation 



3A 3/2 = A(B + -> 7i+K°) + V2A(B+ -> vr ^) 

« v/2 ^ f -f \A{B + - tt+tt )! e^/ 2(( 5 EW _ e *r) , (1) 

where A3/2 is an isospin amplitude parametrizing the A/ = 1 transition B — ► (7rif) 7=3/2, 
e *«P3/2 is a strong-interaction phase, and e n is the weak phase associated with the quark 
decay b — > uus. The second relation is strictly valid in the SU(3) flavor-symmetry limit; 
however, the factor f K J f n = 1.22 ± 0.01 accounts for the leading (i.e., factorizable) 
corrections to that limit. The crucial new ingredient in (|I|) with respect to the corre- 
sponding relation used in Ref. || is the presence of the parameter <5 EW accounting for 
the contributions of electroweak penguin operators. We have recently shown that in the 
SU(3) limit this parameter is real (i.e., it does not carry a non-trivial strong-interaction 
phase) and calculable in terms of Wilson coefficients and electroweak parameters [[Tj| . 
The result is 

S EW = (1 - k) — ^ = 0.63 ± 0.11 , (2) 

where a = 1/129 is the electromagnetic coupling at the weak scale, A = 0.22 is the 
Wolfenstein parameter, Rf, = \~ 1 \V u b/V c b\ ~ 0.41 ± 0.07 p2[ , and k ~ 0.05 accounts for 



factorizable SU(3)-breaking corrections. The derivation of this result uses the fact that 
the relevant electroweak penguin operators are Fierz-equivalent to the usual current- 
current operators Q\ and Q 2 of the effective weak Hamiltonian for B —>■ ttK decays [ |1~5|] , 
and that in the SU(3) limit the isospin amplitude A 3 / 2 receives a contribution only from 
the combination (Qi + Q 2 ), but not from the difference {Q\ — Q 2 )- 

As in the original GRL method, we must rely on the dynamical assumptions that 
\A(B+ -> tt+tt )! = \A(B- -> 7r-7r°)| and A(B+ -> n + K°) = A(B~ -> n-R ). Whereas 
the first relation follows from the fact that only the current-current operators contribute 
to B ± — > tt^tt decays (electroweak penguin contributions can be neglected in this case 
|23|), the second one assumes that there are only negligible contributions proportional 



to the weak phase e* 7 to the amplitude for the decay B + — > tt + K°, which thus can 
be taken to have the simple form A{B + -> tt + K°) = e in e i<l,p \A{B + -> ir + K°)\, where 
e m is the weak phase of the leading top- and charm-penguin amplitudes, and e^ p is a 
strong-interaction phase. Possible contributions to this amplitude proportional to the 
weak phase e n are indeed expected to be very small, because they could only come from 
up-quark penguins or annihilation topologies [24]. However, this intuitive argument 
could be invalidated if soft final-state rescattering effects were very important [14] , ITS] , 
16], 0, 0, gfj. We stress, therefore, that the assumption A(B + -> rr + K°) = A(B~ -> 



7T~K°) is a working hypothesis of our method, which must be tested independently. A 
necessary condition for the validity of this assumption is the absence of a sizable direct CP 
asymmetry in the decays B ± — ► tt ± K°. If we write A(B + — ► n + K°) oc e L ^ p {e m +e %1 e ir, e a ), 
where e a < 1 measures the strength of possible rescattering contributions and e %v is a 
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strong-interaction phase, then acp ~ 2e a sin 7 sin?/. Since the global analysis of the 
unitarity triangle prefers values of 7 such that sin 7 = 0(1), and since sin?] is unlikely 
to be small because without sizable strong phases there would not be a rescattering 
contribution in the first place, a small experimental value for the asymmetry would be 
a strong indication that our working hypothesis is justified. 
Let us define the amplitude ratios 



£3/2 



V us f K V2\A{B^ 



7T 7T 



Wir \A(B+^7T+K°)\ 

V2\A(B ± ^tt ^)! 



\A(B+ -> vr+if )! 



(3) 



which under the assumptions stated above can be determined experimentally through 
time-independent rate measurements via 



£3/2 



V2 



K 



Br(5 H 



7T 1 7T 



Br(£" 



7T 7T 



ON 1 1/2 



Br(5+ - 
Br(fi ± 



tt+K ) + Br(5- 



1/2 



(4) 



A future measurement of r + 7^ r_ would signal direct CP violation in the decays B ± — > 
7r°K . At present, preliminary data reported by the CLEO Collaboration p5| imply 
[k{ r + + r 2 )] x / 2 = 1.46 ± 0.37 and, combined with some theoretical guidance, e 3 / 2 = 
0.24 ±0.06 n Moreover, we define 



'EW 



3 *7 



glz) e 



with z = cos 7, so that 
g(z) 



2z5 E w + 8 2 



ew > 



sin?/' 



sin 7 

5(2) 



(5) 



(6) 



In terms of these quantities, the triangle relation (0) and its CP-conjugate take the form 



l+e 3/2 0(z)e 



i(A0=F^) 



r± e 



(7) 



where A0 = 03/2 — <pp is an unknown strong-interaction phase difference, while the 
phases £± contain both strong and weak contributions. It follows that 



cos(i/j =F 

cos(2?/>) 



£3/2^ (^) 



2e 3 /2^(^) 



= z±(z) 



2(1 



£ 2 (z) 



(8) 
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Combining these results, we find that the allowed solutions for z = cos 7 can be obtained 
from the real zeros of the equation 



16sj /2 g 2 (z) 



X+[Z ) X-(Z, 



(9) 



which, taking into account the z dependence of g(z) and x±(z), correspond to the zeros 
of a fourth-order polynomial in z. 

A simplified analysis can be performed if the phase difference A0 turns out to be 
small or close to 180° - a possibility that can be tested for experimentally. To this end, 
one exploits the following exact relations: 

cos 7 = d E w — 



2e 3/2 (cos A<p + e 3 / 2 5 E w) 

^3/2 " 

The global analysis of the unitarity triangle prefers values of 7 in the range 47° < 7 < 



r 2 + — r 2 _ = sin 7 sin A<p . (10) 



105° which would imply sin 7 > 0.73. Then the second relation can be used to 
obtain a reasonable estimate and upper limit for sin A<p. If it turns out that sin A</> is 
small, corresponding to a situation where |A0| « 0° or 180°, one can set cosA</> = ±1 
in the first relation to obtain 

(1±£ 3 /2<W) 2 -§K + r ~)+4/2 nu 

COS 7~ Or f-M_L g A S ' 

^3/2 (±1 + ^3/2(JEWj 

which determines cos 7 up to a possible two-fold ambiguity. From fllOD , it follows that 
a criterion for the validity of this approximation is that the deviation of cosA0 from 
±1 be less than the uncertainty in the product £3/2<5ew> i-e. min(|A0|, \A(f> — ir\) < 



2A(£ 3 / 2 ( i >ew)- With present uncertainties on the parameters e 3 / 2 and Jew 5 which are 
unlikely to be improved much in the near future, this implies min(|A0|, \A<f> — 180°|) < 
17°. With the current experimental values for the various parameters, and in the absence 
of independent experimental results for r + and r_, the relations (10) do not yet provide 



for a useful estimate of cos 7; however, they may become valuable with more precise 
measurements. It is remarkable that even in the case r + = r_, i.e., in the absence of 
direct CP violation in B ± — > ■n°K ± decays, cos 7 can be determined using relation (|TT|), 
which becomes exact in that limit. 

In practise, the determination of 7 using @) or (|l0|) is limited by experimental as 
well as theoretical uncertainties in the extraction of the parameters r±, £3/2, and Jew- 
Let us illustrate the situation with a realistic example. Assume that the true values of 
the parameters are 7 = 76° (the center of the region preferred by the global analysis), 
£3/2 = 0.24 and Jew — 0.63 (the current central values), and that the strong phase 
difference takes the value A0 = 20°. It then follows that r + rs 1.18 and r_ m 1.04. Let 
us assume that we can measure the values of these parameters with some errors given by 
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Figure 1: Real solutions for z = cos 7 obtained from (Q) in a simulation of 
1000 experiments with Gaussian errors as specified in the text. The correct 
value is z = cos 76° « 0.242. 



AE3/2 = 0.04, A5ew = 0.09, and Ar± = 0.05. We do not anticipate that £ 3 / 2 and 5ew 
will soon be known with an accuracy much better than today, because these quantities 
are affected by theoretical uncertainties such as the estimate of SU(3)-breaking effects. 
We thus assign a 15% error to them [Z7||. The assumed error on the amplitude ratios 
r± corresponds to a measurement of the corresponding ratios of branching ratios with 
a precision of about 10%. In this example, the approximate value for cos 7 obtained by 
setting cos A0 = 1 in (0) is cos 7 « 0.26 ± 0.14(r±) ± 0.09(5 EW ) ± 0.09(e 3/a ), which is 
close to the correct value cos 7 ~ 0.242. We have quoted the various sources of errors 
separately. It is apparent that the precision in the measurements of the ratios r± is the 
limiting factor of our method. The approximate solution obtained with cos A0 = — 1 is 
COS7 ~ 1.13±0.19(r±) ±0.12(5 EW ) ±0.10(63/2), which is excluded by the global analysis 
of the unitarity triangle. In Figure , we show the distribution of the exact real solutions 
of equation (Q) for 1000 random choices of the Gaussian errors in the various input 
quantities. The solutions where cos 7 ~ ±1 can again be excluded based on the global 
analysis of the unitarity triangle. From the central peak, we obtain cos 7 = 0.24 ± 0.18, 
implying at one standard deviation I7I = (76 ± 11)°. 

To conclude, we have shown that the weak phase 7 = Arg(V^ b ) can be determined 
using time-independent measurements of the branching ratios for the decays B + — > tt°K + 
and B~ — > tt°K~ , as well as of the CP-averaged branching ratios for the decays B ± — > 
tt^K and B ± — > 7r ± 7r°. The new development that makes this method practical is the 
observation that the strong phases of the I = | electroweak penguin and tree amplitudes 
are related to one another by Fierz identities and SU(3) flavor symmetry. SU(3)-breaking 
corrections can be accounted for in the factorization approximation. On the other hand, 
like many earlier proposals our method relies on the dynamical assumption that final- 
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state rescatterings do not induce a sizable contribution proportional to the weak phase 
e n in the amplitude for the process B + — > -k + K° . The validity of this assumption can 
be tested for experimentally by searching for direct CP violation in this decay. 
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